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Along with such widely known forms of thermo -chemical process- 
ing of alloys as cementation, nitridation, and cyanidation, certain new 
processes have acquired wide application in recent years. These 
include such processes of diffusion saturation of the surface of alloys 
as saturation with metals (aluminizing, chromizing, molybdenumizing, 
etc. ) and saturation with metalloids (boronizing and siliconizing). 

The extensive industrial use being made of these processes is due 
to the high grade of physico-chemical properties acquired by products 
as a result of diffusion saturation of their surfaces by metals, metal- 
loids, or their compounds. Diffusion saturation of the surface of alloys 
with metals and metalloids is not only an effective means but, in many 
cases, the only possible means for solving such modern technological 
problems as the protection of the surface of alloys from oxidation at 
high temperatures, increasing their durability and resistance to corro- 
sion, providing special electrical and magnetic properties on the sur- 
face, etc. 

Despite the availability of a fairly extensive literature devoted to 
the study of individual questions of diffusion saturation of the surfaces 
of alloys with metals and metalloids (see, for instance, the Literature 
Cited 1,2,3,4,s ), much less attention has been paid to calculating the dif- 
fusion constants in chemico -thermal processing. 

The most thorough investigations in this area are referenced 1,6,7,n 
in the Literature Cited'*'. At the same time, a study of the laws of the 
diffusion of various elements when they are used to saturate the surface 
of metals and alloys makes it possible not only to broaden our under- 
standing of the nature of phase transformations taking place in alloys 
during their processing, but also makes it possible to compute correctly 
the production rate of the diffusion saturation of the surface of alloys 
by means of metals and metalloids. Such computations are of special 
interest in connection with the increased use, in production, of vapor, 
gas, and liquid methods of saturating the surface of alloys by metals 
and metalloids since they make it possible to intensify the saturation 
process. 

The diffusion constants for alpha-phase chromium used for the 
chromizing of iron in a gas medium have already been computed by 


*Given the framework of the present article, we are not citing here 
the extensive literature available on general questions of diffusion. 

This can be had, for instance, from the bibliography cited by Gertsriken 
and Dekhtyar 8 . 
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Hoar and Groom 9 while the corresponding figures for a liquid medium 
are given by Campbell et al 10 . However, the authors of these papers 
in using the basic diffusion equation 

i = Zvy/Dt 

(where £ is the thickness of the diffusion layer; v, a constant which 
depends on concentration and temperature; t, the duration of the dif- 
fusion) conditionally accepted the value of v as being equal to 1. The 
computed values of D and Q then become approximately as follows: 

D^ 00 ° = 14.7 • 1 0~ 10 cm 2 / sec, 

Q = 57 Kcal/g*at (Hoar and Groom 9 ), 

D^ 00 ° = 25. 5 * 1 O” 10 cm 2 / sec, 

Q = 58 Kcal/g-at (Campbell et al 10 ). 

Buckle 11 , using a method he had worked out for computing the dif- 
fusion coefficient for a multiphase layer formed on the surface of a 
sample, utilized the data given by Galmiche 12 on the diffusion-chroming 
of iron in a gas medium, and computed the diffusion constant for the 
alpha- and gamma-phase of chromium, taking into account the concen- 
tration factor of the parameter of v. 

dJj 00 ° = 50* 1 0~ 10 cm 2 / sec, 

D“ 00 ° = 0. 92 • 1 O -10 cm 2 / sec. 

The present article examines the method for determining the dif- 
fusion coefficient in cases when a multiphase diffusion layer is formed 
on the surface zone of the metal. 

The calculation of the coefficient of diffusion was made by using as 
example the chromium -saturation of iron in a gas phase. 

Figure 1 presents schematically the distribution of concentration 
in a sample when, in the process of diffusion, several different phases 
are formed, k = 1,2, . . . n 4 1. 

In determining the diffusion coefficients for various phases in 
the form of a concentration curve, we shall assume that the coefficients 
are constants within the limits of their respective phase and that the 
concentration on the outer surface of sample z 0 is not dependent on time. 
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Figure 1. Schematic Distribution of Concentration 
of Diffusion in a Heterophase Medium 


Concentrations zk-i and ck, on the phase boundaries, are also assumed 
to be known and not to be a function of time. These values can be deter- 
mined in a diagram giving the state of the respective system if one 
assumes that the shift of phase boundaries occurs quasi-statically. 
Generally, the values of ck, zk, and values of the jumps or discontinu- 
ities Ak = ck - zk are determined experimentally. Thus, the concen- 
tration c k (x, t) within phase k satisfies the following equations: 

Bk c xx( x » t) = c t" ( x » *•) 

Yk-i^) < x < Yk(t) 

c k [yk-i(t).t] = zk-i 

c k [yk(t)»t] = c k 

for which the following expression is the solution 

c^ 1 (x, t) = Akerf (x/2>/Dkt) + B k (2) 

where 



A k “ 




erf (b k -3./VDk) " erf(bk/\/Dk) 
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Bk 


ckerf (bk-i/y/Dk) - zk-i erf (bk/\/5k) 
erf (bk-i/y/Dk) - erfCbk/v'Dk) 


6 k = z k-i - ck • 


Within the limitations outlined above and in the presence of natural 
initial conditions c 11 "*" 1 (x, 0) = 0, the boundaries separating the differ- 
ent phases shift— as we know— during the process of diffusion, according 
to the parabolic law. 

It is, therefore, assumed in equations (1) and (Z) that 

y k (t) = 2b k \/F . (3) 

The index k in equations (1) to (3) assumes the value of 1,2, ... n + 1, 
with yo(t) - 0 (b 0 = 0); while when k = n + 1 , we must regard 
yn+i(t) = 00 (t>n + i = °°)> c n + i = 0 (see Figure 1). The mass balance 
on an arbitrary boundary, moving according to the law of x = y(t), can 
be expressed by equation 

dm/dt = - Dc x |y (t), tj - c[y (t), tjdy (t)/dt (4) 

where dm/dt is the resulting flow across the phase boundary in question 
in the direction of an increase in the value of the coordinate x. The 
first and second terms on the right side of the above equation express, 
respectively, the pure diffusion flow and the flow resulting from the 
crossing of the phase boundary by the given substance because of the 
boundary shift. The quantity of matter m^, which at a moment in time 
t is found beyond the border which separates phases k and k + 1, moving 
in accordance with the law of x = yk(t) = 2bkVtT, can be computed by 
the following two methods (numerically, the value of mk equals the 
shaded area in Figure 1): 
t 

m k = - Dk Jcx[yk(t), tldt - ckyk(t) 

o t 

= - Dk + i Jcx + 1 fyk(t), tl dt - Z k yk(t) • (5) 

o L J 

The equivalence of the two methods is the result of the continuity of the 
flow across the phase borders. 


Using equations (2) and (3), we can then compute the following: 

m k ./~D k a, exp( -bk/Dk) . (M 

Zv^V » k erf (b k /VDk) - erfftk^/VDk) ' k k ’ 
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m k _ j£ k+i s exp (-bk/Dk+ 1 ) 

Zy/t f tt 1 erf (bk + i/VDk+ 1 ) - erf(bk/ V^k+i) Zk k ' 


By replacing 
mk 


O ' i ^ • 

:-i/2 \/F + zk-ibk-i _ / bk - bk-A 
m k /2y/T + ckbk ex Py Dk J 


which, solved for Dk, becomes 

H-i + z k-lYk-i\ 
mk + ckyk / 


4Dkt - (yk - Yk-i)/ ln ^- 


(8) 


(9) 


This expression can be used for determining the coefficient of diffusion 
in various phases. To find the coefficient Dk, it is sufficient to know 
the position of a given phase boundary at any given point in time t, and 
the relationship of the integral diffusion flows across these boundaries 
as of that particular moment, 0k = ( m k + l + z k-i Yk- i)/ ( m k + c kYk)- 


According to Figure 1, the value of 9k represents the relationship 
between the areas of sections of the diffusion zone lying below the a^^ 
line and the ak line, respectively. It can be calculated directly if the 
distribution concentration in the given sample has been determined 
exp e r im e ntally . 

Formula (9), which for an area adjacent to the surface of the sam- 
ple (k = 1) takes the form of 


4Dit 


yf/ln 


( 


m 0 

m : + cjyi 



( 10 ) 


can be utilized for determining the coefficient of diffusion in all zones 
except the last one, since when k = n + 1 (i. e. , when yn + 1 00 and 

c n + 1 — 0), expression (9) becomes identical. Formula (7) can be used 
to determine the value of D n + : . With k = n, it gives 


V^e^erfce = ^ , 

m n + z nYn 

P = Yu 

2VDn+ it ' 


(ID 


The right side of this equation can also be determined graphically, if 
the course of the concentration curve in the last zone is known. The 
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function y/n x exp (x 2 ) erfcx is a monotone function, the value of which 
increases and tends to approach 1 for increasing values of x. 

Following is a table of the values of this function, which can be 
used to find the required value of £ in equation (11): 


X 

• xe x2 arfcx 

X 

y/n • xe x2 erfcx 

X 

>fn • xe x2 erfcx 

X 

y/n • xe x2 erfcx 

0. 1 

0. 159 

0. 6 

0. 603 

1.2 

0.804 

2.2 

0. 917 

0. 2 

0.286 

0.7 

0. 652 

1.4 

0.839 

2.4 

0. 929 

0. 3 

0. 390 

0.8 

0. 693 

1.6 

0.866 

2.6 

0. 936 

0.4 

0.475 

0.9 

0. 720 

1.8 

0.887 

2.8 

0. 944 

0. 5 

0. 545 

1. 0 

0. 757 

2. 0 

0. 904 

3. 0 

0. 950 


Expressions (6) and (7) can serve as a basis for another method of 
determining the diffusion coefficient of a multiphase medium. This 
method— unlike the one described above— does not require knowledge of 
the exact progress of concentration in the sample under investigation. 
Since expressions (6) and (7) are equal, we obtain the following recur- 
rent equation 

5, 2 V Dk + it . exp ( -yk/ 4 Dk+ it) 

Kt 1 Yk ‘ erf (yk+ 1 /2 V D k + it) - erf (yk/ 2 V^k + it) 

Ak , (12) 

which makes it possible to determine the value of Dk + i if Dk h as already 
been determined and if the position of the phase -boundaries at the time- 
instant t—as well as the concentration level at these points— is known. 

Thus, in order to find all the coefficients, it is sufficient— besides estab- 
lishing the value of S^, Ak, and yk, already mentioned— to know just the 
value of Dj. 

To find Di, convenient use can be made of the formula which deter- 
mines the total mass m 0 of the substance which had diffused its way into 
the body of the sample. We must bear in mind, however, that in a num- 
ber of cases it is possible to determine m 0 without defining the form of 
the concentration curve (e. g. , by means of weighing the sample if sat- 
uration took place in a gas or liquid phase, or by measuring the total 
y -radioactivity of the sample with the use of radioactive isotopes). 

The equation for m 0 derived from formula (7) for the case when 
k = 0 has this form: 


= 6 k 2 ^ /D k t _ exp ( - ykMPkt) 

yk ’ erf (yk/2V^>kt) - erf (yk - i /2 V^kt) 
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( 13 ) 


By using formula (12), the coefficients D 2 , D 3 , . . . D n + i can be 
computed by the graphic method if we first draw up a family of functions 
exp ( - x 2 )/x (erfObc - erfx) to determine the values of x and OL at the 
intervals required (x = yk/2 • V^k + it, OL = yk+i/yk ^ 1 for the left 
half of the equation, and x = yk/2 V^kU OL = yk~i/yk ^ 1 for the right). 

By using the first of the methods described above, we were able, 
with the aid of formulas (9) to (11), to determine the coefficients for 
the diffusion of chromium in iron (with a 0. 03% C content), for OL and 
y phases, saturating the iron with chromium in the gas phase 13 . 

Figure 2, below, shows the results of determining the concentra- 
tion of chromium by layers in the diffusion zone at temperatures of 
900°, 1000°, 1100°, and 1200°C— the duration of the process being 3 
and 6 hours. 



Figure 2. Changes in the Concentration of Chromium in the Surface 
Zones of Iron Which Underwent a Chromium Saturation Process 
in a Gas Medium at Temperatures from 900° to 1200°C for a 
Period of Three Hours (1) and Six Hours (2) 


Chromium concentration by layers was determined by the spectral 
method 14 . The nodal points of the concentration curves (identified by 
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arrows in Figure 2) correspond to the boundaries of y-*Ci phase- 
transition at an isothermic point in the process. The upper part of 
the curve corresponds to chromium distribution in the Qf-phase, while 
the lower part corresponds to the y-phase (at annealing temperature). 
Chromium concentration on the boundary of the Q? and y phases, obtained 
as a result of these experiments, corresponds satisfactorily to the 
results shown in the diagram of the Fe = Cr system state 15 . 

The results of computing the chromium diffusion coefficient in the 
Of and y phases are shown in Figure 3, below. 



Figure 3. The Temperature Factor in the Coefficient 
of Diffusion (D, cm 2 /sec) of Chromium in the &- 
State of Iron (Curve 1) and the y-State (Curve 2) 


Processing the experimental data makes it possible to present 
equations for the chromium diffusion coefficient in the Qt and y phases 
in iron as follows: 



(14) 

(15) 


Received for publication: 7 Oct 1963. 
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